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ABSTRACT 

This paper reviews the status of singular value loop-shaping as 
a design paradigm for multivariable feedback systems. It shows 
that this paradigm is an effective design tool whenever the prob- 
lem specifications are ’spacially round’. The tool can be arbitrarily 
conservative, however, when they are not. This happens because 
singular value conditions for robust performance are not tight 
(necessary and sufficient) and can severely overstate actual 
requirements. An alternate paradigm is discussed which over- 
comes these limitations. The alternative includes a more general 
problem formulation, a new matrix function (. i , and tight conditions 
for both robust stability and robust performance. The state of the 
art currently permits analysis of feedback systems within this new 
paradigm. Synthesis remains a subject of research. 
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1. INTRODUCTION 

Ever since the basic work of Nyquist [l], Bode [2] and others, the classical 
approach to feedback design has followed a frequency domain perspective. We 
are given a plant described by a rational transfer function, C(s), and wish to 
design a rational compensator, K(s), such that the closed loop feedback system 
is stable and meets certain performance and robustness requirements. 

As is well known, the stability requirement imposes structural constraints 
on certain transfer functions of the closed loop system, e.g. a Nyquist encircle- 
ment count for the function det(I+GK) [3]. Likewise, the performance require- 
ment imposes magnitude constraints on certain other transfer functions. In 
particular, for disturbances and commands reflected to the output loop- 
breaking point, the (output) sensitivity function 

5(s)^[/ + £?( S )A'(s)]-» (1) 

must be small for all frequencies, s=ju, where the disturbances and/or refer- 
ence commands are large. 

The third feedback design requirement — tolerance for uncertainty -- also 
imposes magnitude constraints on transfer functions. A common requirement 
in this case is that the "complementary (output) sensitivity function" 

T(s) = G(s)K(s)[I+G(s) K(s))~' (2) 

must be small for all frequencies where so-called unstructured multiplicative 
model uncertainties are large [4]. 

For classical single-input single-output (SISO) systems, the meanings of 
"small" and "large” in these statements are, of course, understood in terms of 
the absolute values of the respective complex-valued functions at each fre- 
quency. Hence, SISO designers working in the frequency domain have viewed 
the design problem as one of shaping the (Bode) magnitude plots of sensitivity 
and complementary sensitivity functions to be small enough to meet design 
specifications. Indeed, because |S(s)| 1/ | GK(s ) | whenever iGA'(s)! is large, 

and |r(s){ « |G/r(s)| whenever | GK(s)\ is small, the shapes of these magni- 
tude functions are intimately tied to the shape of the loop transfer function 
GK(s), and the entire design process is often referred to simply as "loop- 
shaping". 

Over the last few years, the loop-shaping process has been successfully gen- 
eralized to multi-input multi-output (MIMO) design problems [4]. Key 
ingredients of the generalization include -- 1) the use of singular values as 
appropriate measures of magnitude for matrix-valued transfer functions, 2) the 
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development of formal mathematical conditions -which guarantee stability 
robustness and performance robustness of MIMO feedback systems in terms of 
these magnitude measures, and 3) the development of certain modifications of 
existing design procedures (e.g. LQG/LTR) which help to synthesize desired mul- 
tivariable loop shapes. 

Design experience with these new results shows that loop-shaping is an 
effective MIMO design paradigm for problems whose specifications can be 
reduced to "spacially round" requirements on S(s) and T(s) alone. Unfor- 
tunately, many design issues which arise in MIMO problems cannot be usefully 
expressed in this form. This paper describes some of these latter design issues 
and develops a design framework and certain recently developed tools which 
promise to deal with them more effectively. 

The paper begins in Section 2 with a brief review of the MIMO loop-shaping 
process. It then discusses some of the issues which are not easily handled by 
this process in Section 3. Finally, a more general design framework and associ- 
ated research results are presented in Section 4 which promise to address these 
issues more effectively. Section 5 provides concluding comments. The paper 
presents no basic new theoretical results and should be viewed only as a brief 
look at the current status of this branch of frequency domain MIMO design. 
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2. THE MIMO LOOP-SHAPING GENERALIZATION 

Our generic multivariable feedback design problem is illustrated in Figure 
1. The loop consists of a plant and a compensator in a unity feedback arrange- 
ment. The plant can be any element from a set of plants characterized by a 
nominal operator, G. and a perturbation operator, dG. These operators are 
modelled by rational transfer function matrices G(s ) and 6G(s), respectively. 
Similarly, the compensator, K is modelled by K(s). The design problem is to 
find a A"(s) which makes the feedback loop internally stable for all possible 
plants G(s)+6G(s) and causes it to respond well to external signals such as 
commands, r(s), disturbances, d(s), and sensor noise, n(s).* 


2.1. A Formal Design Problem 

This MIMO design problem can be formalized by specifying precise 
mathematical statements for the qualitative performance objectives above, and 
by specifying a set of external signals, r(s), d(s) and n(s). and a set of plants, 
G(s)+6G(s), over which these objectives must be achieved. We will start with 
some very simple specifications for these three elements of the design problem. 
More complex situations are treated later. 


2.1.1. Performance Objectives 

As formal performance objectives, we will require that the error signals 
from the feedback loop, e(s) = y(s)-r(s), be sufficiently small in the Lg-norm 
sense. That is, w r e require 

lH 

/ e(t) T e(t)dt < 6 (3) 


u £ 
i ! S — 


for all external signals and all plants in sets yet to be defined. The scalar 6 sets 
the desired level of performance. 


2.1.2. External Signals 

We also use the Lg-norm to characterize the set of external signals. For 
the time being, this set will consist of a single vector-valued signal, r(s), with 
d(s) s 0, and n(s) = 0. The signal r(s) will consists of all time functions which 
can be generated by passing the functions 

v = with Ihll* * i (a) 

through linear systems with specified frequency responses, i.e. 

r(s) = w r (s)r](s) (5) 


Note that these signals belong to a unit ball in Lg w r hich has been distorted in 
frequency content by the system w r I. with transfer function iv r (s)I, to 
represent the spectral content of commands. Without loss of generality, ir r I is 
assumed to be stable and to have a stable inverse. 


2.1.3. Magnitudes of Transfer Matrices 

The above choices of performance objectives and signal sets imply that 
singular values are the appropriate measure of magnitude for matrix transfer 
functions. This follows from the operator norm induced by Lg-functions [5]. To 

* Throughout the paper, x(t) end x(s) will he used to designate time functions end their La- 
place transforms, respectively. 
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illustrate, suppose that <5G is identically zero. Then 

e(s) = S(s)r(s) (6) 

where S(s) is the feedback system's nominal sensitivity matrix defined in equa- 
tion (l). The performance objective for this case requires that 

sup |! e || 2 = sup || Sw r I 77 1 | 2 
a ** 

= lisw r i;i 2 , 2 

= sup a[S(;£j)^ r (;'u)3 <6. (7) 

U 

Here the function d[] denotes the largest singular value of its matrix-valued 
argument. It follows that the performance objective is achieved (nominally) if 
and only if the Bode plot of this singular value for S(ju) lies entirely below the 
Bode plot of 6/ |u) r (jf cj) j over the entire frequency range 0 ^ u < «*> . 

2. 1.4. Set of Plants 

The operator norm induced by external L 2 -signals also provides a con- 
venient (and reasonable [6] ) way to express the set of plants. As a starting 
point, we will consider G+6G to be generated by unstructured multiplicative per- 
turbations at the output, i.e. 

[C?(s)+dC(s)j = [/ + I(s)]G(s) , (B) 

where L(s) is a perturbation matrix "which satisfies 

L(s) = v l (s)A(s) (9) 

for a specified stable invertible w^I and some arbitrary stable* operator A with 
induced norm less than or equal to unity, i.e. 

I! A ;i 2-2 = sup a[A(?'y)] £1 (10) 

U 

Note that this characterization makes L norm-bounded in the sense that L 2 - 
signals in the unit ball, as shaped in frequency content by Wx -1 I, produce out- 
puts whose L 2 -norms are less than unity. 

2.2. Robust Stability and Robust Performance 

Given these specifications for performance objectives, external signal sets, 
and plant sets, w r e now seek a compensator K(s) which satisfies two require- 
ments. First, it must achieve stability for all elements in the defined set of 
plants. This property will be called "robust stability". Second, it must satisfy 
the stated performance objective for all signals in the defined signal set and for 
all plants in the plant set. This latter property will be called "robust perfor- 
mance". The process of finding such a compensator is greatly facilitated by two 
important analysis results. The first of these is the following theorem, proven in 
[4], which provides necessary and sufficient conditions for robust stability. 

Theorem 1 

Suppose that the nominal feedback system in Figure 1 is stable (i.e. it is 
stable when 6G(s) s 0). Then the perturbed system is stable for all <5£?(s) 
defined by equations (9)-(l0) if and only if 


• The stability requirement on A can be relaxed to the condition that G end (1+L)G have 

the same number of unstable modes. 
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o[‘w L {jd) T{ju)] < 1 for all u, (ll) 

where T(s) is the complementary sensitivity function defined by equation 

( 2 ). 

This result follows from the fact that the feedback loop remains stable if and 
only if the function det [I+(I+L)GK] = det[f+GK]det[I+&WLT] remains non- 
zero along the j cj-axis (and therefore in the right half plane) for all A. 

The second important analysis result provides a sufficient condition and a 
separate necessary condition for robust performance: 

Theorem 2 

Suppose that the feedback system in Figure 1 is robustly stable. Then the 
perturbed system satisfies performance objective (3) for all signals defined 
by (4)-(5) and all 6G(s) defined by (B)-(10) if 

a[50'y)^r0'u)] < <5 [ 1 -o[w L (ju) T(ju)] ] for all u , (12) 

and it satisfies the objective for all signals and all plants only if 

a[S( 7 'y)uy( 7 'y)] < <5 [ (7«) T(ju) ] for ail u , (-13) 

where jt[] denotes the smallest singular value of its argument. 

These conditions follow directly from the perturbed sensitivity function 

5(s) = / + [/+I(s)]G(s)tf(s) 

= S(s)[H-L(s)T(s)]-' (14) 

which shows that performance is maintained in the face of L(s) whenever the 
nominal performance requirement. a[Su; r ]^ 6 from equation (7), is tightened 
sufficiently to offset the amplification of the factor [I+LT ] -1 . Equations (12) 
and (13) simply reflect the worst and best case values this amplification can 
take. Note that these equations reduce to a single necessary and sufficient con- 
dition ■whenever the complementary sensitivity function is "spacially round”, i.e. 
when its condition number satisfies 

k[T ] = o[T]/2.[T] * 1. (15) 

Even -without this property, however, the sufficient condition (12) alone is not 
unduely conservative. This is so because the factor 1 —olyuiT] must be positive 
for robust stability and is typically designed to be 0.5 or greater to provide some 
design margin. On the other hand, the factor is never greater than 

unity. Hence, the true robust performance requirement is typically overstated 
by less than a factor of two. 

The significance of Theorems 1 and 2 is that they define robust stability and 
robust performance solely in terms of acceptable magnitudes for the nominal 
functions a[ T(j y) ] and a[ S(j u) ]. In particular, we saw from equation (11) that 
robust stability is achieved if and only if 

° [T(ju)] < \w L (ju) ~ f0rallU - 

Similarly, equation (12) shows that a linear combination of d[ F ] and a[S] must 
be small enough to assure robust performance, namely 
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a[S(ju) 3 + 1 o[ T{j’J) ] < -r — f-ry for all o . (17) 

\vj r (ju)\ L |l^rO") I 

As in the SISO case, therefore, the design problem defined by the specified per- 
formance objectives, external signal sets, and plant sets in Sections 1.1. 1-1. 1.4 
reduces to one of shaping Bode plots of sensitivity and complementary sensi- 
tivity magnitudes such that they lie below specified bounds over the entire fre- 
quency range. 

2.3. Synthesis Methods 

While our simple MIMO design problem is thus fundamentally the same as a 
SISO problem, the actual process of shaping MIMO functions is more difficult 
than shaping SISO ones. This difficulty has been demonstated over the years by 
a variety of attempts to generalize SISO synthesis concepts (e.g. inverse Nyquist 
methods applied to diagonally dominant systems [7], direct-Nyquist and Bode 
methods applied to characteristic loci [B], root-locus methods applied to mul- 
tivariable functions [9], etc). 

It turns out that some of the most effective methods of shaping MIMO loops 
use modern optimization-based synthesis tools. For example, modified versions 
of the LQG problem can be used effectively to achieve trade-offs between singu- 
lar values of S(s) and T(s) across frequency [4,10]. More recent H«-methods 
have also been developed to synthesize compensators which directly minimize 
sup a [S(jcu)], or sup o[T(j w)], or singular values of weighted augmented combi- 

tt U 

nations supu[5(ju) j w(jv)T(j '->)] [11,12,13]. However, no methods have been 

U 

developed to date which directly optimize sums of weighted singular values, as 
might be suggested by equation (17). 

These various synthesis methods will not be discussed further here. Rather, 
it is our objective to identify design issues which do not lend themselves readily 
to the MIMO loop-shaping generalization even if the associated synthesis prob- 
lems could be solved effectively. 
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3. MORE COMPLEX PROBLEMS 

The design problem described above can be made more general and more 
useful for design purposes by including more complex performance require- 
ments, signal sets, and/or plant sets. An obvious addition, for example, is to 
include disturbances d(s) in the external signals. If these are generated as in 
equation (5) by an L 2 -ball shaped in frequency content by w*I, then a robust per- 
formance result analogous to equation (17) applies, with |u;<i(s)| replacing 
|u; r (s)j. If both r(s) and d(s) exist simultaneously, then i u* ! + j ui r | replaces 
\w r j alone. 

Another obvious addition is to include sensor noise, n(s). Let this signal 
again be generated by an L 2 -ball shaped by w n I. With other signals zero, this 
input drives the error e (s ) via 

e(s) = (I+L)GK[I + (I+L)GK]~ l vj n r)(s) 

= (I+L)T[I + LT]-'w n r,(s). (IB) 


which (after some algebra) leads to the following sufficient condition for perfor- 
mance over all n(s) and all L(s): 


o[T(ju)] < 


6 | vji (J u) j + j «*»(>«)! (l+|t«iO'“)!) 


for all u 


(19) 


Note that this constraint on afF] is tighter than the stability robustness con- 
straint (16). Hence, it can be used in place of (16) to cover both design require- 
ments. 


3.1. Difficulties 

While the above two signal additions can be accomodated quite easily in the 
MIMO loop-shaping setup, there are other generalizations which cannot. In fact, 
any performance requirement, signal set, or plant set which is not spacially 
round as seen at the feedback system’s output can cause difficulties.* To illus- 
trate this, we will consider the following more general specifications: 

Performance Requirements 

Filtered versions of e(s) must be Lg-norm bounded for all external inputs 
and all plants, i.e. 

|| c i! a < <5 with c(s) = B«(s)e(s), (20) 

■where ^(s ) is the transfer matrix of a stable invertible operator, W t , which 
shapes the spacial directions as well as the frequency content of e (s). 

External Signals 
r(s) generated by 

r(s) = W r (s)rj(s ) , (21) 

where matrix RV( S ) shapes the spacial directions and frequency content of 
the commands. 


• Technically, any specification of performance requirements, signal sets, and plant sets 
which are either all spacially round at the output or all spacially round at the input can be 
accomodated. In the latter case, the design would be done with S and T defined at the in- 
put. 



- 8 - 


Set of Plants 

Unstructured multiplicative perturbations at the output defined by 

L(s) = W Lo (s)h(s)W Li (s) , (22) 

where ^(s) and ^(s) shape the spacial directions and frequency content 
of inputs and outputs of the unit-norm perturbation A, respectively. 

M1M0 loop-shaping results for these specifications can be developed in a manner 
completely analogous to Section 2. Under the assumption that the nominal 
feedback system is stable, the principal results are the following: 

(1) Nominal performance is satisfied if and only if 

a[ S(j'S) W r (Ju)] < 6 for all u (23) 

(2) Stability is robust if and only if 

o[Wn{ju)T{j :j) H’io (; tS) ] < 1 for all u (24) 

(3) Performance is robust if (not only if) 

5[ W*SW r ] < y [ 1 -a[ W Li TW L 0 ] J for all a. (25) 

The result which causes difficulties here is (25). It is derived via the following 
manipulations: 

i{s) = JJ'g (s)S(s)W T (s) r)(s ) 

= W a S[I + LT]-'W T -n(s) 

= W’gSWV Wr-'il+LTY'Wr’ois) 

= (^5FV r )(^ 1 lK r )- 1 [/ + A(^ i T^)]-U^ 1 WV)r?(s) (26) 

The last of these equations shows that the tightest singular value bound which 
we can place on the transfer matrix from r)(s) to c(s) for all L(s) is 
a[fV, / ( 1 -g\_Wl\. TWlo])- This bound leads directly to Condition 
(25), and thus, (25) is the weakest sufficient condition for robust stability which 
can be established via singular values. 

Unfortunately, (25) can be arbitrarily conservative — requiring much 
tighter nominal performance than is actually necessary to assure robust perfor- 
mance. Examples which illustrate this difficulty abound [6]. Indeed, (25) is 
often found to be impossibly tight in light of stability robustness constraints. 
The potential for this conservatism exists whenever the condition number 
K [^Lo^r] is large.* This occurs when the external signals and/or the plant per- 
turbations are not spacially round. 

3.2. Examples 

To conclude this section, we briefly consider two examples which give rise 
to problem specifications which are not spacially round. 

Example 1: Perturbations at the Input 

Suppose that the set of plants. G+5G, is described not by (8) in Section 2 
but rather by 

• To a lesser degree, the potential also exists when /c[ Wn T Wig ]) is large. As in Section 2, 
however, the factor 1— aX Wn TVt'i 0 ] is typically designed to be 0.5 or greater, end thus, its 
contribution to conservatism is not excessive. 
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lG(s) + 6G(s)] = G{s)U + L(s)] (2?) 

rvi 

with L(s) satisfying assumptions (9)-(l0). Then a few manipulations show 
that the corresponding multiplicative perturbation at the output is given by 

L(s) = G(s)L(s) G(s)-' (2B) 

Thus, L(s) satisfies Specification (22) with 

V u (s) = /(s)C.(s)-» and W u (s) = G.(s) w L (s) / f (s) (29) 

where G.(s) is a minimum phase stable version of G(s) (with all unstable 
poles and transmission zeros reflected about the jcj- axis), and / (s) is a low 
pass filter which makes G» _1 proper at high frequencies. It follows that this 
specification of plant perturbations, which is spacially round at the input, 
will not be round at the output if the condition number of G{ju) is large. 

Example 2: Disturbances at the Input 

Suppose that the external signals consist of disturbances du( s ) entering at 
the input and generated by an Lz-ball shaped with w u I. These disturbances 
can be represented by the following equivalent* disturbances reflected to 
the output: 

d(s) = G.(s)d u (s) (30) 

This disturbance corresponds to Specification (21) with D’V(s) = i u u (s)G»(s). 
Again, this external signal is not round whenever the condition number of 
G(j n) is large. 

Both of these examples are, of course, common occurances in MIMO feedback 
design, and there are many others. 


•Equivalent in the sense that closed loop responses y (s) will have the same magnitudes for 
all s=ju. 
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4. AN ALTERNATE MIMO DESIGN FRAMEWORK 

The loop-shaping limitations described above can be overcome with an 
alternate design framework which has been developed in the last few years 
[6,14,15]. This alternate framework consists of a new problem description, a 
new measure of magnitude for matrix transfer functions, and certain key 
analysis and synthesis results. These various elements are described briefly in 
this Section. 

4.1. Problem Description 

The new problem description is illustrated in Figure 2. It consists of a very 
general "plant", P, whose outputs and inputs comprise three pairs of vector 
variables. The first pair of variables consists of measured outputs, y(s), to be 
used for feedback and control inputs, u (s), to be commanded by the compen- 
sator. K. The second pair consists of performance variables, c(s). and external 
input signals, T)(s). Finally, the third pair of variables consists of input signals, 
ft(s), to the perturbation operator, A, and outputs from this operator, £ 0 (s), 
which feed back into the plant. The design problem is to find a compensator 
which keeps the signals s(s) in the unit L z -ball for all signals rj(s) in the unit- 
ball and all stable unit-norm operators A. 

This problem description is very general because the internal structure of 
P can be chosen to represent many different problem specifications. One 
example of this internal structure is shown in Figure 3. This figure corresponds 
to the problem specifications in Section 3.1. P is seen to include the usual 
input-output description of the real plant G. but it also includes the weighting 
operators W„ and W r , which shape performance variables and external signals, 
as well as and Wr* which shape the plant perturbation. The types of exter- 
nal signals (whether r(s), d(s), n(s) and/or others) are also defined by the 
internal structure, as are the locations of perturbations (whether at outputs, 
inputs, and/or elsewhere). Various examples of internal structures of P for 
other problem specifications can be found in [6], 

4.2. Analysis Results 

Beyond mere generality, Figure 2 is important because it comes equipped 
with a non-conservative necessary and sufficient condition for robust perfor- 
mance. In order to describe this new condition, w-e first close the compensator 
feedback loop in Figure 2 to get the closed loop system in Figure 4. The opera- 
tor F(P,K) in this figure has a 2x2 block-structured transfer function matrix 
F(s) whose blocks are defined in terms of the original 3x3 partition of P(s) as 
follows: 

Fij(s) = Fy (s) + P i3 (s)[/ -K(s)P 33 (s)]-' K(s)P 3j (s) i.j = 1.2 (31) 

Suppose that this system is stable. Then the following results apply: 

(1) Nominal performance is satisfied if and only if 

(7 [ Fzzij «j) ] < 1 for all co (32) 

(2) Stability is robust if and only if 

a [.F n (ju)] < 1 for all co (33) 

(3) Performance and stability are robust if and only if 

fx.[F(jv)] < 1 for all co 


(34) 
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•where /i[] is a function to be defined shortly. 

The first two of these results are self-evident. Result(32) follows by definition — 
performance is satisfied for A(s) s 0 iff the induced norm ||F 22 !| 2 „ 2 is less than 
or equal to unity. Result(33) follows from the stability condition with the A-loop 
closed -- namely det [/-AF n ] must remain non-zero on the j oaxis for all A. 

The more significant result is (34). This follows again from the definition 
that performance is robust if and only if the t/r} transfer matrix -with the A-loop 
closed remains stable and norm-bounded by unity, i.e. iff (33) is satisfied and 

a[.F 22 + Fz\ (I -&Fu)~ l &Fiz] < 1 for all u and all A (35) 

Notice that this last norm-bound is also a necessary and sufficient condition for 
continued stability if we chose to connect a second norm-bounded perturbation, 
say A 0 (s), across the t and tj terminals of Figure 4 (to see this, compare the 
form of (35) with our other stability conditions (ll), (24), and (33) ). It follows, 
therefore, that robust performance is equivalent to robust stability in the face 
of two perturbations, A and A 0 , connected around the system F(s) in the diago- 
nally structured arrangement shown in Figure 5. 

These observations bring us to the function ji[]. This function was defined 
in [16] to provide a magnitude measure for the smallest block-structured per- 
turbation which will make a system unstable. The full definition of y for com- 
plex matrices is the following: 

I det [ I—5XM ] = 0 

min <5 | for some X = cfiap[Ai,A 2 , • • • ,A m ] 
i with a [ Aj ] < 1 for all i 

In words, this equation defines /u[] to be the reciprocal of the smallest value of 
scalar <5 which makes the matrix I-6XM singular for some A’ in a block-diagonal 
perturbation set. Notice that this definition reduces to the conventional singular 
value in the absence of structure (i.e. when the number of blocks, m, in X is 
one). For this reason, y has been called the "structured singular value". Note 
also that the value of y depends on the number of blocks in the structure as 
well as on the dimensions of these blocks. Technically, therefore, y's argu- 
ments should include not only matrix M but also a multi-index which describes 
the structure. By convention and for sake of notational simplicity, this latter 
dependence is suppressed. 

It is clear from this definition that y can be applied to the transfer 
matrices of Figure 5 to test whether de£[/-diaj(A, A 0 )F] remains non-zero 
along the jio- axis. This establishes tight conditions for robust stability with 
respect to the two perturbation blocks, and equivalently, tight conditions for 
robust stability and performance (Condition (34)). Formal details of this argu- 
ment are given in [14]. Notice, however, that the definition of y is not limited to 
the 2x2 diagonal perturbation block structure in Figure 5. It can be used to 
test stability with respect to any number of diagonal blocks. This makes it pos- 
sible to establish robust stability with respect to plant sets which are charac- 
terized by several unstructured perturbations, and simultaneously, to establish 
robust performance with respect to several performance requirements*. 
Indeed, the only limitation on perturbation structures testable via y appears to 
be that each perturbation block must be allowed to be complex-valued. 
Research to remove this remaining potential source of conservatism is under 
way [17]. 

• The meaning of several 'performance blocks' is not entirely clear. For such cases, y< 1, 
implies that the i-th performance requirement is robust with respect to the reel perturbs- 
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4.3. Numerics for /* 

Like singular values, ji is useful for practical numerical analyses as well as 
for theoretical ones. First-generation computer algorithms have been 
developed to evaluate the function for fixed complex matrices. When used 
repeatedly, these algorithms can generate Bode plots of /t over frequency for 
matrix transfer functions such as F(s). This provides a practical Bode-like 
analysis test of stability/performance robustness for any given candidate 
design. 

To date, /^-algorithms are based on the the following inequalities, proven in 

[16]: 

• max j\[ UM ] | sJ /*[/</] rmn a[D M D''] (37) 

where 

U = diag [ U\, Uz. • ■ • , U m ] Ui unitary 

D - dia.g [dily, dzlz, ■ ■ • .dnln] rij scalar 

Reference [16] shows that the left hand side of inequality (37) is tight and thus 
prorides a potential way to compute fi. Unfortunately, the implied maximiza- 
tion over the block-structured unitary matrices U has many local maxima. The 
right hand side of (37) is also tight, at least for structures with three or few r er 
blocks. Its implied minimization over the block-structured scaling matrices D 
is convex and thus provides a much nicer problem for numerial search solu- 
tions. For this reason, current /i-algorithms are based on the right had side of 
(37). The issues posed by four or more blocks in the structure remain under 
study, and further improvements in algorithms are forthcoming [19]. 

4.4. /i-Synthesis 

Progress has also been made in the development of formal synthesis 
methods for the design framework in Figure 2. These methods seek to design 
compensators, K to stabilize the nominal system, P, and to minimize jx[F]. 
While complete solutions of this problem are not yet available, an iterative 
scheme has been invented which yields useful answers f 19]. The iterative 
scheme exploits the fact exhibited in (37) that fx is a scaled version of d with 
block-structured scaling matrix D. This fact suggests the following iteration: 

Step(l) 

Fix an initial estimate for D{jtS) 

Step(2) 

Solve a a-synthesis problem to find a stabilizing K(s) which minimizes 
sup d[D(J u)F(J 'S)D(j u)~ l ] 

U 

Step(3) 

Evaluate the /i-properties of this solution against Condition (34), and 
Step(4) 

Stop if Condition (34) is satisfied. Else replace the D(j’J) estimate with the 
D(jo) values computed as part of the /^-calculation in Step(3) and return 
to Step(2) 


lions A as well as the This is a stronger requirement than actually necessary. 
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These iterations are practical primarily because the cr-problem in Step(2) has a 
numerically tractable solution. This a-solution has been completed only 
recently [15], and is itself a significant step forward in MIMO design. It encom- 
passes the various special //.-problems from Section 2 as special cases and 
provides state-space-based computational algorithms which can handle design 
problems of significant engineering size. Detailed descriptions of the solution 
are left to [15] . 
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5. SUMMARY AND CONCLUSIONS 

This paper has provided a brief status review of the singular value loop- 
shaping design paradigm for multivarible feedback control systems. It has 
shown that this paradigm is useful for design problems whose specifications of 
external signal sets, plant sets, and performance requirements are spacially 
round at the plant output or. by duality, spacially round at the plant input. For 
such problems, it is possible to write tight necessary and sufficient analysis con- 
ditions for robust stability and also reasonably tight conditions for robust per- 
formance. Both conditions take the form of Bode-like magnitude bounds on 
singular values of sensititivity and complementary sensitivity matrices of the 
feedback system. 

Unfortunately, when problem specifications are not spacially round, the 
singular value conditions for robust performance can be arbitrarily conserva- 
tive, leading to highly overspecified design requirements. Design problems in 
this category abound. An alternate design paradigm is discussed which over- 
comes these limitations. The alternate paradigm includes of a more general 
problem formulation, a new matrix function /i , and tight necessary and sufficient 
conditions for both robust stability and robust performance in terms of this new 
function. Under current state of the art, fj. can be calculated numerically and 
thus provides an effective analysis tool for existing candidate designs. However, 
direct synthesis techniques to design compensators w'hich stabilize and also 
minimize (J. are not yet available. Research to provide this capability is under- 
way. 
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